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UBIQUITOUS TORSIONAL MOTIONS IN TYPE II SPICULES 
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ABSTRACT 

Spicules are long, thin, highly dynamic features that jut out ubiquitously from the solar limb. They 
dominate the interface between the chromosphere and corona and may provide significant mass and 
energy to the corona. We use high-quality observations with the Swedish 1-m Solar Telescope to 
establish that so-called type II spicules are characterized by the simultaneous action of three different 
types of motion: (i) field-aligned flows of order 50-100 kms -1 , (ii) swaying motions of order 15-20 
kms -1 , and (Hi), torsional motions of order 25-30 kms -1 . The first two modes have been studied 
in detail before, but not the torsional motions. Our analysis of many near-limb and off-limb spectra 
and narrow-band images using multiple spectral lines yields strong evidence that most, if not all, 
type-II spicules undergo large torsional modulation and that these motions, like spicule swaying, 
represent Alfvenic waves propagating outward at several hundred kms -1 . The combined action of the 
different motions explains the similar morphology of spicule bushes in the outer red and blue wings of 
chromospheric lines, and needs to be taken into account when interpreting Doppler motions to derive 
estimates for field-aligned flows in spicules and determining the Alfvenic wave energy in the solar 
atmosphere. Our results also suggest that large torsional motion is an ingredient in the production of 
type-II spicules and that spicules play an important role in the transport of helicity through the solar 
atmosphere. 

Subject headings: Sun: chromosphere - Sun: corona - Sun: oscillations - Sun: magnetic fields 



1. INTRODUCTION 

Between the photosphere and corona lies the chromo- 
sphere, a region of relatively cool plasma that is most 
conspicuous in the hydrogen Ha Balmer line. The up- 
per chromosphere is dominated by spicules, thin jets of 
chromospheric plasma that reach heights of 10 000 km 
or more above the photosphere. Although spicules 
were described already by Secchi in 1877, understand- 
ing their physical natur e has progresse d only slowly (re- 
views bv IBeckersI fl968: Sterling 2000). The launch of 
the Hinode satellite (|Kosugi et al.l l2Q0"7h and the com- 
bined use of adaptive optics and image post-processing 
(|van Noort et al.l [20051 ) in ground-based observing have 
revolutionized our view of spicules. 

The re are (at least) two types (|De Pontieu et all 
I2007H ). Type I spicules reach heights of 2-9 Mm, 
show up- and downward velocities of 10-30 kms -1 , 
and have lifetimes of 3-10 minutes. They probably 
correspond to on-disk dynamic fibrils caused by shock 
waves that propagate upward into the chromosphere 
(lHansteen et al. 2006: lRouppe van der Voort et al. 2007; 
iDe Pontieu et al.ll2007af ) (see also[Suematsu et al. 1 9951 ). 
Type II spicules reach larger heights at velocities of or- 
der 50-100 kms -1 . They have shorter lifetimes, of 
order 100s, and usually only exhibit upward motion 
before their rapid fading in the chromospheric lines 
in which they are detected. On the disk they ap- 
pear as rapidly moving absorpt ion features in the blue 
wings of chromospheric lines (E angange n et al.l 2008; 
iRouppe van der Voort et all 120091 ). In this paper we fo- 
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cus on this intriguing class of features. 

Our previous studies suggest that type II spicules 
represent impulsively accelerated chromospheric ma- 
terial that is continuously heated while it rises 
(|De Pontieu et al.l I2QQ9L 1201 ll) . The cause of the heat- 
ing and acceleration are unknown, but a magnetic pro- 
cess such as reconnect ion and /or flux emergence is 
most likely (IDe Pontieuetal.l2QQ7bl : ISterling etaLll201tt 
IMartmez-Svkora et al.ll2011f ). 

Type II spicules show other motions in addition to 
radial outflow. In the Ca II H line they are seen to 
sway to-and-from transversely with amplitudes of order 
10-20 kms -1 a nd periodicities of 100- 500 s, suggesting 
Alfvenic waves ([De Pontieu et al"ll2007cf ). The continu- 
ation of these motions in transition region and coronal 
lines suggests that th ey may help drive the solar wind 
([Mcintosh et al.ll201lD . 

Other types of m otion are less well established. 
iSuematsu et al.l (|2008f ) reported that some spicules ap- 
pear as double t hread s with evidence of s p inning motion. 
ICurdt fc Tianl ([201lh and ICurdt et al.l (j2012l ) suggest 
that similar spinning explains the tilts of ultraviolet lines 
in so-called explosive events producing larger-scale macro 
spicules. Spectral-line tilts were noted earlier in observa- 
tions at the lim b and also attribute d to spicule rotation 
(IBeckersI 1 19721 ). or not interpreted ([Houtgast fc Namibal 
I1979D . 

In this paper we definitely confirm the indications for 
twisting spicular motions by unequivocally detecting tor- 
sional spicule modulation in limb spectroscopy and imag- 
ing spectroscopy with unprecedented spectral and spatial 
resolution. 
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Figure 1. Examples of large torsional motion in off-limb spicules in CallH. Top row shows the spatial variation of the CallH profile, 
with the right hand panel of each pair a slit-jaw image with the off-limb intensities enhanced with a radial filter. These show spicules 
protruding upwards from the limb which is to the right. The lefthand panels are corresponding CallH spectrograms. The black line in 
the slit-jaw images shows the slit location. The spectrograms are dominated by linear features that are often tilted from the horizontal 
dispersion direction, a clear signature of torsional motions, typically over 20-40 kms -1 . Selected spicules are marked by black ticks in the 
slit-jaw images and white ticks in the spectrograms. Bottom row shows examples of large swaying motion in off-limb spicules in CallH. 
Each panel shows the temporal evolution of the CallH profile. Despite significant superposition of multiple oscillatory signals there are 
often episodes of quasi-sinusoidal oscillation. For guidance we have overplotted sinus curves with the amplitude ("amp" in kms -1 ) and 
period (P in s) specified in each legend. 



2. OBSERVATIONS 

We analyze various data sets obt ained with the 
Swedis h 1-m Solar Telescope (SST, iScharmer et al.l 
2003a) on La Palma: slit spectroscopy with the TRI- 
Port Polarimetric Echelle- Littrow spectrograph (TRIP- 
PEL, Kiselm an et al.l 120111 ) and imaging spectroscopy 
with the CRisp Ima ging SpectroPolarimeter (CRISP, 
IScharmer et al.l 12008). All observations used adaptive 
optics. 

With the TRIPPEL spectrograph we observed the 
CallH line at 396.8 nm with nominal spectral resolu- 
tion of 1.6 pm. A slit-jaw camera of the same type was 
slaved to the spectrum camera. The data were binned to 
0.49 pm spectral and 0.068" spatial pixel sizes. The spec- 
trograms were corrected for dark curren t, gain variations, 
and spectrograph distortions following Langangen et al. 
(2007]). We use a 55min duration, 0.91s cadence se- 
quence of CallH spectrograms taken on 2009 October 8 
in excellent seeing. The slit was set parallel to the limb 
at various heights above the limb ranging from 3" to 10". 



The exposure time was 800 ms. 

CRISP contains a dual Fabry-Perot interferometer and 
allows for fast (<50 ms) wavelength tuning within a spec- 
tral region. High spatial resolution and precise alignment 
between the sequentially taken images for different tun- 
ings is achieved with t he i mage restoration technique of 
Ivan Noort et al.l (|2005[ ) and iRouppe van der Voort et al.l 
(|2009f ). In section [3] we present results from two CRISP 
data sets: a 17min sequence at the ultra-high cadence 
of 0.44 s sampling Ha (656.3 nm) only at ±1204 mA 
(±55 kms -1 ), registering the limb near AR11230 on 2011 
June 11, and an earlier 36min, 17 s cadence sequence of 
finely spaced Ha and Ca II 854.2 nm profile samplings 
registering the limb on 2010 June 27. 

3. RESULTS 

Our various observations all show clear evidence of 
torsional motions on very small scales, of order 0.7" or 
less. We first demonstrate these in the TRIPPEL spec- 
tra. Since the slit was oriented approximately parallel 
to the limb, it crossed many spicules. In Figure Q] (top 




Torsional motions in spicules 



Observations 



Simulation, best fit 



Observations 



" 20 




Simulation, no swaying Simulation, no torsion 




-100 



-50 50 
Velocity [km/s] 



1-OC00 



-50 50 
Velocity [km/s] 



100 



Figure 2. Comparison between observations and Monte Carlo 
simulations in X—y displays. The best-fit solution has field-aligned 
velocities of 60 kms -1 , torsional velocities of 30 kms -1 and sway- 
ing motions of 15 kms -1 with periods of 100-300 s. Significant 
discrepancies occur for larger swaying, more torsion, no swaying, 
and no torsion. 

row) we show two X — y cuts, selecting moments when 
the slit was located at least several arcseconds above the 
limb to avoid the enormous line-of-sight superposition 
at and just above the limb. The slit-jaw images indeed 
show many distinct spicules that seem at least partially 
resolved. Each spectral X — y panel contains multiple 
cases in which a spicule appears as a thin bright streak 
across the spectrum, with a small tilt from the horizontal 
wavelength direction and a substantial offset from nom- 
inal line center. We detected such behavior throughout 
the 55-min sequence at all sampled heights. The tilted- 
streak morphology indicates relative redshift on one side 
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Figure 3. Comparison between observations and Monte Carlo 
simulations in A— t displays. The best-fit solution has field-aligned 
velocities of 60 kms -1 , torsional velocities of 30 kms -1 and sway- 
ing motions of 15 kms -1 with periods of 100-300 s. Significant 
discrepancies occur for wave periods that are much higher or lower 
than these values. 

of a spicule, blueshift on the other. This reversal in trans- 
verse motion is the signature of torsional spicule motion. 
In addition, the substantial offsets of the tilted streaks 
from nominal line center can be understood as the su- 
perposition of swaying motion and a projection of field- 
aligned flow onto the line-of-sight. 

Next, we examine how our data can be used to sepa- 
rate and constrain the magnitudes and other properties 
of these three types of spicular motion. Figure |T] (bottom 
row) shows the temporal variation of the CallH profile. 
Although significant line-of-sight superposition occurs, 
we find clear evidence of quasi-sinusoidal episodes in the 
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Figure 4. Example of sign change in spicule Doppler tilt in the ultra-high cadence (0.44 s) ±55 kms -1 CRISP Ha data. Lefthand panel 
shows an Ha wing image at maximum spicule length. Center panels contain sequences of four Doppler images at different sampling times 
(in seconds, blueshift bright). The righthand panel show space-time diagrams measured along the spicule axis. Dashed guidelines illustrate 
a propagation speed of 285 kms -1 . 



Doppler shifts. To get an idea of the typical amplitudes 
and periods of such oscillations we overplot simple func- 
tions v = vo + a sm(2ir(t—to) / P) and find that amplitudes 
of order 20-40 kms -1 and periods of order 100-200 s are 
common. Swaying and torsional motions can both pro- 
duce such sinusoidal signature in A — t plots, but they 
have a different appearance in a time sequence of succes- 
sive A — y plots, with swaying motions showing wiggling 
of the entire streak-like feature whereas torsional mod- 
ulation shows up as a tilted streak that shrinks into a 
more vertical feature and then becomes tilted with oppo- 
site sign. Detailed inspection of many such plots showed 
that temporal undulation as in Figure [T] (bottom row) is 
sometimes associated with swaying and sometimes with 
torsional modulation. 

Our data inspection thus reveals evidence for three dis- 
tinct types of spicular motion. Given the limitations of 
our data sets (in particular the enormous line-of-sight su- 
perposition at the limb and variations in seeing quality), 
we find that the best method to derive statistical proper- 
ties of these different motions is to compare the observa- 
tions with Monte Carlo si mulations. We adapt met hod- 
ology used prev i ously by iDe Pontieu et al.l ((200 7c) and 
iMcIntosh et al.l ([2008): we consider TV spicules and im- 
pose on each a field- aligned flow v p , a torsional motion 
Vt, and a swaying motion v s . Inspired by our earlier 
analyses of type II spicules we assume that, during their 
lifetime T, they continuously grow with constant veloc- 
ity v p along a straight path, inclined by an angle a from 
the vertical and with an azimuth angle /3, until they fade 
rapidly from view. The torsional and swaying motions 
are assumed to have periods P t and P s , respectively, with 
random phase. 

Figures [2] and [3] shows examples of comparisons be- 
tween observed A — y and A — t plots and results from 
these Monte Carlo simulations for different parameter 
choices. We made many such plots with large parame- 
ter variation and defined a best-fit solution by selecting 
parameter combinations that reproduce the observed be- 
havior the best. Per spicule this best-fit choice takes v p 
randomly from a Gaussian distribution around 60 kms -1 
with standard deviation a = 10 kms -1 (similar to what 
is found by Pereira et al., in prep.), T from a Gaussian 
distribution around 120s with <j = 30s (Pereira et al., in 



prep.), a from a Gaussian distribution around 20 degrees 
with cr = 10 degrees (based on slit-jaw images), f3 from a 
uniform distribution over 0-360 degrees, v t from a Gaus- 
sian distribution around 30 kms -1 with a = 10 kms -1 , 
v s from a Gaussian distribution around 15 kms -1 with 
<T = 5 kms -1 , and both P s and Pt from a uniform distri- 
bution over 100-300 s. Comparison of the A — y and A — t 
panels for this best-fit solution with the observations in 
the top left panels of Figures [2] and [3] and in Figure [T] 
shows that the best-fit solution reproduces, statistically, 
the appearance of both types of data, with a similar mul- 
titude of slightly tilted streak-like features in X — y and 
similar sinusoidal swings in X — t. 

How well defined are the best-fit parameters, given 
the large number of free parameters in our simulations? 
To answer this question we ran multiple simulations in 
which we changed the distribution for only one parame- 
ter, keeping all others fixed. Examples are given in the 
lower rows of Figures [2] and 03 In the first panel of the sec- 
ond row (Fig. [2j) the mean swaying amplitude is doubled 
and gives an overall zig-zag pattern with too many ex- 
tremes. Similarly, doubling the mean torsional amplitude 
gives streaks that are too wide in velocity. Removing the 
swaying (bottom-left panel) produces too little zig-zag 
motion in X — y. No torsional motion gives streaks that 
are too narrow in velocity. Similarly, in the X—t panels of 
Fig. [3j doubling the frequencies to periods of 50-100 s of 
the swaying or torsional motion leads to too many wig- 
gles, whereas lower frequencies (periods 300-600 s) yield 
too few wiggles. Note that it is more difficult to de- 
termine the perio ds since these often exce ed the spicule 
lifetimes (see also lDe Pontieu et al.ll2007cl h 

In summary, our Monte- Carlo analysis provides reason- 
ably well-defined constraints. In order to reproduce the 
appearance of our limb spectra the torsional and sway- 
ing motions should be of order 30 and 15 kms -1 , respec- 
tively, with periodicities of order 100 — 300s. 

The next issue is whether the transverse swaying 
and torsional modulations represent propagating waves. 
Since our limb spectra sample only one height above the 
limb at a time, we address this question with CRISP 
imaging spectroscopy in H<x The ultra-high cadence se- 
quence permits the construction of Doppler images by 
subtracting the images taken in the red and blue wings 
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(corresponding to ±55 kms -1 Dopplershift). In these we 
often observe very fast propagation of the Doppler sig- 
nal. An example is shown in Figure |4] where the apparent 
phase speed is 285 kms -1 . Values in the range 100-300 
kms -1 are common. This is the order of magnitude ex- 
pected for the Alfven speed in stru ctures with densities 
of order 10 10 cm -3 (Beckers 196 81) and magnetic field 
strengths of order 10-30 Gauss (jCenteno et al .1120101 ). 

Determining whether these Dopplershift modulations 
represent swaying or torsional motion is not straightfor- 
ward. In our Dopplergrams, torsional modulation will 
show up as a black-and-white pattern across a spicule 
only when the sum of the swaying motion and projection 
of the field-aligned flow happens to be zero. The large 
offsets from line center in Figure [T] suggest that such can- 
cellation does not occur often. Direct separation of the 
torsional and swaying modes is impeded by the combi- 
nation of very high phase speed, sparse wavelength sam- 
pling, and relatively low cadence, i.e., lack of simultaneity 
between the red- and blue- wing samplings. Nevertheless, 
the black-and-white pattern of the spicule at t = 34.32 s 
in Figure H] seems direct evidence of significant torsional 
motion. Its propagation at about 300 kms -1 provides 
further support that the observed rotational motions are 
a signature of torsional Alfven waves that propagate out- 
ward along spicules. 

4. DISCUSSION 

Our observations and analysis lead to the conclusion 
that spicules undergo three different types of motions at 
the same time: field-aligned flow aligned with the spicule 
axis, swaying that moves the spicule as a whole trans- 
verse to its axis, and torsional motion around its axis. 
The superposition of these motions complicates the inter- 
pretation of spicule spectroscopy and imaging consider- 
ably. Our Monte Carlo simulations served to disentangle 
them and indicate that our observations are compatible 
with field-aligned flows of order 50-100 kms -1 , swaying 
of order 15-20 kms -1 , and torsional motion of order 25- 
30 kms -1 . We also found evidence that the swaying 
and torsional motions are both signatures of Alfvenic 
waves with periods of order 100-300 s with propagation 
along the spicule axis at phase speeds as high as 100- 
300 kms -1 . Our data are inadequate to decide whether 
only upward propagation occurs, or whether propagation 
in both directions or even standing waves are as impor- 
tant for the torsi onal motions as they seem t o be for 
swaying motions (Okam oto fc De PontieulfeOllI ). 

Our results impact several long-standing issues. First, 
any interpretation of observed spicular motion is likely 
wrong if not all three kinetic modes are properly ac- 
counted for. The classic reports of field-aligned outflows 
of or der 20-30 km s -1 from Doppler shifts at the limb 
(e.g., Beckers 1972f ) did not include torsional and swaying 
contributions. This introduces additional uncertainty in 
the historically reported values for field-aligned spicular 
flows, which are already plagued by poor spatio-temporal 
resolution. 

Second, juxtaposition of the three modes of mo- 
tion explains the appearance of spicular features in 
the blue and red wings of Call 854.2 nm and Ha 
close to the limb. Such images show rows ("bushes") 
of near-vertical absorption features that are compara- 
ble to the on-disk "rapid blue excursions" (RBEs) of 



iRouppe van der Voort et all (|2009l ). When these are 
compared between near-simultaneous red and blue outer- 
wing samplings, the red-shifted and blue-shifted features 
appear with remarkably similar, though not identical, 
morphology and orientation. Precise fine-scale overlap 
is rare but the bush patterns often appear roughly the 
same. An Ha example pair is shown in the lower panels 
of Figure [5] Such outer- wing similarity has been inter- 
preted as presence of similar field-aligned flows point- 
ing towards and away from the observer ([Becker si fl 9 72), 
but the asymmetry of near-limb line-of-sight projection 
should then cause large morphological inequality. It is 
much more likely that large- amplitude transverse mo- 
tions dominate the visibility close to the limb of outer- 
wing features and that they cause morphological equal- 
ity through near-cospatial and alternating red and blue 
Dopplershift modulation. The upper panels of Figure [5] 
indeed show spicular Doppler tilts in the central spec- 
trum panels, much like the off-limb ones in Figure [TJ 
This suggests that on-disk RBEs also undergo torsional 
motion. 

Our observ a tions confirm the speculation of 
iCurdt fc Tianl (|2011[ ) that torsional motions likely 
occur on smaller scales than their SUMER observations 
of macrospi cules and explosive events (or the so-called 
swirl s, see IWedemever-Bohm fc Rouppe van der Voort 
2009). Our results also sup port the early wo r k by Avery 
(1970) and suggestions by Pasac hoff et al.l (|1968f ) that 
upward propagating Alfven waves explain indications 
of high-speed propagation of Doppler signals along 
spicules. The presence and/or observed periods of prop- 
agating torsional Alfven waves can in principle provide 
information about the thermodynamic a nd magnetic 
properties of the guid ing structures (e.g., iRouth et al.l 
2010: IVerth et al.ll2010f ). Such idealized models are not 
yet a good representation of actual, highly dynamic 
spicules but may illustrate the potential diagnostic value 
of spicular wave observations. 

The ubiquity of torsional spicule motion also pro- 
vides support for scenarios of spicule formation in which 
non-linear coupling of torsional Alfven waves to other 
wave modes on expanding flux tubes leads to sig- 
nificant field-aligned flows that drive the plasma up- 
wards ([H ollweg et al.l 119821: Sterling & Hollwej] 119841: 
| De Pontieu fc Ha erendel 19981: iKudoh fc Shibatal 119991 : 
iMatsumoto fc Shibatal 120101 ) . However, none of these 
models can explain both the large observed field-aligned 
flows and the rapid heating to at least transition-region 
temperatures. This suggests that while torsional mo- 
tions may play a role in providing upward momentum in 
spicules, other components of spicule formation are still 
missing. 

While our Monte Carlo simulation assumes the pres- 
ence of only the m = torsional mode, the observa- 
tions show also evidence for more complex torsional 
modes with m > (see example at y=34" in top left 
panel of Figure [T]). Our results thus may provide sup- 
port for the recent hypothesis (Ivan B allegoo i3en et al.l 
1201 ll : lAsgari-Targhi fc van Ballegooijenll2012f ) that com- 
plex torsional motions are generated in the interaction 
between convective motions and photospheric flux tubes 
(see, e.g., t he ubiqitous vor ticy reported in recent sim- 
ulations by iMoll et all l2011f ) , and that they may be re- 
sponsible for heating of the chromosphere and corona. 



6 



De Pontieu et al. 





f 1 « 




-2 -1 1 -20 -10 10 20 
x [arcsec] Doppler [km/s] 



-2-10 1 
x [arcsec] 



-2 -1 1 -20 -10 10 20 
x [arcsec] Doppler [km/s] 



-2-10 1 
x [arcsec] 



S 40 



" 35 




10 



20 30 
x [arcsec] 



40 



50 



10 



20 30 
x [arcsec] 



40 



50 



Figure 5. Two examples of CRISP Call 854.2 nm spectral behavior close to the limb which is to the left (in top row). The slit-like 
sampling crossed a bush of near-limb spicules as shown by the black line in the Ha context images in the lower row. In the top row, each 
Call A — y diagram is accompanied by two Call slit-jaw images that have the location of the slit marked by a black line, at left a wing 
image, at right a Doppler image. The Call spectrograms and images were taken 17s apart, the righthand trio first. The Ha images were 
taken within a few seconds of each other 34 s later. 



The observed torsional motions imply that spicules, 
given their ubiquity, play a major role in the trans- 
port of helicity through the solar atmosphere. In addi- 
tion, our results suggest that the energy flux carried by 
Alfvenic motions into the corona and heliosphere may 
be twice as large as the previ ous estimates that were 
based on swaying mo tions only (jDe Pontieu et al.ll2QQ7d : 
iMcIntosh etaIll2Qllh . 
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